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ABSTRACT 


A  system  Is  presented  that  combines  the  automated  planning 
and.  optimization  functions  in  machining  processes.   The  planning 
function  is  performed  by  a  systematic  analysis  of  the  stated  require- 
ments of  the  finished  part  in  the  light  of  information  on  available 
machining  facilities  and  raw  materials.  The  optimization  phase 
utilizes  a  mathematical  programming  model  to  take  into  account 
various  costs  and  constraints  under  alternative  machining  conditions. 
A  gradient  or  'hill-climbing'  algorithm  is  shown  to  be  a  convenient 
optimization  technique  for  this  class  of  problems.   Implementation  of 
the  system  is  Illustrated  in  some  detail  for  the  specific  case  of 
•the  face  milling  process. 
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INTRODUCTION 


In  conjunction  with  the  recent  advances  in  computer 
technology,  modern  industry  is  rapidly  moving  toward  a  cor.pletcly 
automated  factory.  There  has  been  considerable  progress  in  the 
areas  of  computer-aided  design  as  well  as  conputer  controlled 
manufacture.   The  field  between  the  design  and  manufacturing  stages 
'that  Is  receiving  increased  research  attention  is  the  area  of 
automated  planning  and  optimization  of  manufacturing  processes. 
Automated  Manufacturing-Planning  (AMP)  is  a  relatively  new  concept 
and  most  of  its  development  nas  taken  place  in  the  sixties.  In  AMP 
design  data  are  converted  into  manufacturing  instructions  that  can 
"be  used  to  make  the  required  finished  component  on  the  selected 
machine  tool.  The  advantages  of  AMP  are  many.  A  uniformly  high 
level  of  planning  and  optimization  can  be  obtained  regardless  of 
work  load.  This  in  turn  leads  to  shorter  manufacturing  cycles  and 
better  utilization  of  the  available  facilities  and  labor.  Automated 
planning  can  also  incorporate  automated  updating  which  takes  into  con- 
sideration day  to  day  changes  in  available  facilities,  time  standards, 
planning  logic,  etc.  '   • 

The  MILVAP  systen  1 15]  and  the  AUTOPIT  system  Il6j 
are  among  tp.e  first  attempts  for  develop-ent  of  automated  planning 
systems  for  machining  problems.  These  systems  achieve  automated  selectior 
of  speeds,  feeds  and  otlicr  machining  parameters,  selection  of 
suitable  tools  and  also  perform  normal  numerical  control  tape 


preparation  functions  such  as   control  of  tool  paths.   Researchers  at 
IBM  [26]  also  developed  their  ovm  "Automated  Manufacturing  Planning 
System"  built  on  similar  lines.  The  "Regenerative  Shop  Planning" 
proposed  by  Scott  118  ]  may  also  be  considered  a  A2IP  system,  \;ith 
the  qualification  that  it  is -usable  only  when  a  number  of  geometrically 
similar  parts  are  encountered  so  that  the  planning  logic  used  in 
the  manufacture  of  one  part  can  be  preserved  and  reused  in  planning 
of  other  parts. 

Several  German 'Universities  and  Industrial  Organizations  are 
developing  the  EXAPT  system  I19].   It  is  divided  into  three  parts, 
the  first  for  point  to  point  drilling  and  milling,  the  second  for 
straight  cutting  and  contouring  on  lathes  and  the  third  for  milling 
straight  cuts  and  contours.  The  first  two  parts  are  commercially 
available  wliile  the  third  part  is  still  under  development  flA!]. 


The  above  systems  do  not.  In  general,  consider  the  constraints 
on  the  selection  of  machining  parameters  that  may  be  dictated  by  a 
variety  of  physical  considerations  in  the  machining  process.  Moreover, 
these  systems  focus  on  the  planning  phase,  paying  relatively  little 
attention  to  optimization.   They  do  not  incorporate  subsystems 
that  explicitly  optimize  the  various  machining  parameters.   On  the 
other  hand,  optimization  of  machining  parameters  has  been  attacked 
many  times  as  a  separate  problem  by  a  number  of  researchers  including 
Taylor  [20],  Gilbert  [9],  Weill  [25],  Colding  [8],  Brewer  and  Reuda  [5], 
and  Okushima  and  Hitomi  [13].  They  all  utilize  empirical  formulas 
for  expressing  tool  life  as  a  function  of  various  machining  parameters 
such  as  cutting  speed,  feed  and  depth  of  cut  and  most  of  them 
differentiate  these  expressions  to  obtain  the  optimum  values  of  the 
parameters.  However,  there  are  difficulties  in  Incorporating  these 
analytical  procedures  into  an  automated  planning  system.  Machining 
processes  typically  involve  a  large  number  of  variables  which  change 
from  one  job  to  another  and  from  one  stage  of  machining  to  another. 
This  makes  the  purely  analytical  methods,  even  in  their  simplest  form, 
impractical  for  applications  involving  complex  machining  sequences 
Involving  large  numbers of  machines,  tools  and  types  of  machining. 
Besides,  most  of  these  analytical  procedures  ignore  the  physical 
constraints  on  the  process  as  well  as  the  probabilistic  nature  of  tool 
failures. 


The  present  paper  is  an  attempt  at  a  systems  synthesis  of 
the  planning  and  optimization  phases.  The  goal  is  to  arrive  at  a 
physically  Implementable  manufacturing  planning  system  that  is 
'self-optimizing'.   This  idea  is  not  completely  new.   Berra  and 
Barash  [3,4]  have  developed  an  automated  planning  system  for  rough 
turning  operations  and  optimized  it  using  an  empirical  search 
procedure.  A  special  feature  of  their  work  was  the  explicit  consider- 
ation of  the  probabilistic  nature  of  cutting  tool  failures.   In  follow 
on  research,  Batra  and  Barash  [2]  extended  the  work  to  consider  multi- 
tool  set-ups  in  turning  operations.   In  the  present  work,  it  is  shown 
that  several  features  peculiar  to  machining  processes  make  it 
particularly  convenient  to  model  the  optimization  phase  as  a  mathematical 
programming  problem.  This  optimization  procedure  overcomes  most  of  the 
difficulties  present  in  other  analytical  procedures,  and  at  the  same 
time  blends  smoothly  into  the  planning  logic. 


SCOPE  OF  THE  PRESENT  WORK 


A  comprehensive  system  of  automated  manufacturing  planning 
must  Include  several  different  phases  of  planning  and  optimization 
such  as:   selection  of  jobs  to  be  undertaken,  determination  of 
optimum  quantities  of  manufacture,  selection  of  broad  machining 
sequence  and  then  the  detailed  machining  sequence  at  each  stage, 
optimization  of  machining  loading,  routing  and  optimum  combination 
of  job  schedules  for  the  different  jobs  in  the  shop.  In  the  present 
work,  the  main  emphasis  is  on  the  phase  of  planning  and  optimization 
of  machining  sequences.   The  system  suggested  herein  would  thus  be 
subsystem  of  a  more  comprehensive  manufacturing  planning  system. 
However,  in  practical  situations,  this  subsystem  alone  can  be 
extremely  complex.  Moreover,  it  is  this  part  of  a  manufacturing 
planning  system  that  is  most  amenable  to  being  run  as  automated 
' self-optimizing '  system.  There  is,  of  course,  no  ban  on  using  this 
subsystem  by  Itself  in  the  absence  of  a  comprehensive  manufacturing 
planning  system. 


The  system  presented  here  accepts  data  on  the  specifications 
of  the  finished  workpiece,  selects  the  best  raw  material  out  of  the 
available  ones,  determines  the  finishing  and  grinding  allowances 
needed  to  achieve  the  required  specifications,  and  hence  the  overall 
machining  requirements  and  the  exact  amount  of  material  (both  shape 
and  size)  to  be  removed  in  each  stage  of  the  machining.   It  then 
optimizes  the  sequence  of  machining  as  well  as  values  of  the  machining 
parameters  at  each  stage.  A  detailed  description  of  all  machines, 
cutting  tools  and  raw  materials  available  for  machining  is  assumed 
to  be  available. 


•  •       -The  output  of  the  system  consists  of  a  print  out  of  tlie 
details  of  the  selected  raw  material,  machine,  cutter,  selected 
values  of  the  machining  parameters  such  as  speed,  feed,  depth  of 
cut  and  v/idth  of  cut  and  a  print  out  of  the  size  of  the  machined 
piece  after  each  stage  of  machining.   It  also  provides  an  estimate 
of  the  total  time  required  for  the  machining  and  the  total  cost 
of  machining  per  piece.  "    ■  ■ 


A  system  block  diagram  for  the  planning  procedure  is 
given  in  Figure  %.     The  following  abbreviations  have  been  used 
for  denoting  the  various  data  and  output  matrices. 


CDM  Component  Data  Matrix 

CDMFC  Finish  Cut  Component  Data  Matrix 

CDMFG   *  Finish  Grind  Component  Data  Matrix 

CDMRC  Rough  Cut  Component  Data  Matrix 

CDMRG  Rough  Grind  Component  Data  Matrix 

MCM  Machine  Characteristics  Matrix 

MRM  Machine  Rcquire^.ent  Matrix 

RCSIZE  Rough  Cut  Size  Matrix 

RMCM  Raw  Material  Characteristics  Matrix 

TOM  Tool  Characteristics  Matrix  . 


INPUT    Itl   THE    FORM 

OF    DATA  MATRICES 

CaM,   RI-'.CM,   MCM 

AND  TC.M 
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GENERATE  OTHER 
DATA  MATRICES  SUCH  . 
AS  CDHRC,  CDMFC,  / 
CDFiRG,  RCSIZE 
MRM,  ETC. 
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Figure  1.   System  Block  Diagram 
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Figure  1.   System  Block  Diagram   (contd.) 


10 


Select 

best  machine 

and  cutter 


PRINT  ALL  RECOM- 
MENDATIONS, COR- 
RESPONDING, 
COST  &  TIME 


P*-TO  C 


C 


[ 


END 


") 


Figure  *.   System  Block  Diagran  (cont'd) 
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The  input  to  the  systcn  consists  of  the  data  matrices 
CDM,  RMCM,  MCM  and  TCM.  The  CDM  describes  all  the  specifications 
of  the  workpieces  to  be  machined  including  all  dimensions,  tolerances 
and  surface  finish  grades  needed.   In  the  present  systen  it  is 
...-issuned  that  the  machining  of  the  part  can  be  subdivided  into  the 
nllling  of'  a  number  of  plain  surfaces  and  hence  each  row  in  the  CDM 
describes  exactly  one  such  plain  surface.  The  RI-'CM  gives  a  complete 
description  of  the  size,  shape  and  the  material  of  the  blank  (or  raw 
material).   The  MCM  and  TCM  descritie  the  various  characteristics  of  the 
■machines  and  tools  to  be  considered.  The  specifications  given  . 
'  in  MCM  include  the  maaimum  and  minimum  limits  for  spindle  speed 
and  feed  rate  on  the  machine,  horsepower  of  the  driving  motor, 
efficiencies  of  the  transmission  of  the  mach"Jnc  and  the  driving 
motor,  static  stiffness  and  the  cost  rates  for  operators,  setters 
and  overheads.   The  TCM  gives  for  each  cutter,  the  number  of  teeth 
on  the  cutter,  diameter  of  the  cutter,  length  of  the  cutting  edge, 
maximum  force  that  can  be  resisted  by  the  tooth  tip  and  values 
of  average  change  tine  and  average  regrinding  tine  for  each  cutter. 

The  first  matrix  generated  by  the  system  is  the  Machine 
pprtulrer'enf  Matrix.   Thp  surface  finish  de«?lp;natiop  and  th«? 
tolerance  specification  for  the  three  dimensions  of  the  work  surface 
are  considered  jointly  to  determine  the  requirements  of  each 
surface  i.e.,  whether  each  of  the  surfaces  requires  some  or  all  of 
the  operations;  rough  billing,  finish  milling,  rough  grinding  and 
finish  grinding. 


•    -       .  •         .  !  •■ 
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The  decision  on  the  r.achine  rcqulrer.cnts  is  nade  based  on  the 
dimension  which  has  the  tightest  tolerance  specification  or  finest 
surface  finish  designation.  The  Machine  Requirements  Matrix  is 
a  pattern  of  I's  and  O's  such  that  a  '1*  in  the  column  opposite 
a  surface  number  indicates  that  the  particular  machining 
operation  must  be  performed  on  that  surface.  Appropriate  machining 
allowances  are  then  added  to  each  of  the  dimensions  of  the  work 
surface  to  generate  the  matrices  CDMRC,  CDMFC  and  CDMRG  which  describe 
the  siZ3  and  shape  of  the  component  after  the  rough  milling,  finish 

?  milling  and  rough  grinding  operations  respectively.  The  system 
then  checks  if  the  raw  material  size  specified  in  the  RMCM  is 
smaller  or  larger  than  that  specified  by  the  CDMRC.   If  it  is 
smaller,  it  prints  out  an  error  message  and  goes  to  the  next  raw 
material  specification.   If  larger,  the  system  determines  the  amount 
of  material  to  be  removed  from  the  given  raw  material  to  bring  it 
to  the  size  specified  by  CDMRC  and  stores  it  in  the  matrix 
RCSiZE.  The  system  repeats  this  process  for  each  of  the 
acceptable  raw  materials.  Thus,  RCSIZE  is  actually  a  three  dimensional 

•  'array  of  numbers;  each  layer  of  the  array  corresponding  to  one 
raw  material.   The  numbers  in  the  array  are  the  dimensions  of  the 
material  to  be  removed  in  the  milling  operation  for  each  »7ork 
Surf ace. Out  of  all  acceptable  raw  materials,  the  best  is  chosen  by 
the  criterion  of  minimum  total  volume  of  material  to  be  removed. 
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The  system  then  selects  the  work  surfaces  one  at  a  time  and  determines 
the  set  of  'feasible'  machines  and  'feasible'  cutting  tools  for  the 
work  surface  picked.   (By  'feasible'  machines,  we  mean  machines  that 
can  meet  the  machining  requirements  of  the  chosen  work  surface;  the 
definition  of  'feasible'  cutting  tools  is  similar).  Next,  it  chooses 
combinations  of  feasible  machines  and  feasible  cutting  tools  one  by 
one,  determines  the  optimum  cutting  conditions  in  each  case  and 
selects  the  combination  that  yields  the  minimum  optimized  cost.   It 
then  prints  out  the  best  machine,  best  cutter  and  the  optimum  machining 
parameters  for  that  combination  along  with  estimates  of  total  cost  and 
total  time  per  machined  piece. 

The  system  concept  suggested  here  is  quite  general  and  may  be 
used  for  a  wide  variety  of  machining  processes  such  as  turning,  milling, 
shaping,  etc.  Fortunately  too,  the  similarities  in  physical  character- 
istics of  different  machining  processes  are  strong  enough  to  permit 
use  of  the  same  optimization  methodology.  For  instance,  the  functional 
relations  between  dependent  variables  (such  as  tool  life  and  cutting 
forces)  and  independent  variables  (such  as  machining  parameters)  have 
the  same  general  form  for  all  machining  processes.  Specific  formulations 
will  indeed  be  needed  for  each  type  of  machining  and  job  environment, 
but  the  general  operation  and  methodology  of  the  system  will  be  the 
same . 


THE  FACE  MILLING  PROCESS 


Since  the  optimization  phase  requires  consideration  of 
certain  physical  details  of  the  machining  process,  this  phase  will 
be  illustrated  for  the  specific  case  of  a  face  milling  process. 


[-* Lx 


Tigure  ^  Face  Milline 
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For  the  benefit  of  the  unfamiliar  reader,  Figure  "Ogives 
a  very  brief  description  of  the  face  milling  process.  The  face  milling 
cutter  has  its  cutting  teeth  on  the  end  and  is  rotated  by  means  of 
the  spindle  of  the  milling  machine.   Simultaneously,  the  work  table 
moves  In  a  longitudinal  direction,  thus  removing  a  thickness  of 
metal  equal  to  the  depth  of  cut  given.  At  the  end  of  each  longitudinal 
stroke,  a  transverse  'feed'  may  be  given  to  remove  an  adjacent 
longitudinal  strip.  The  machining  parameters  of  this  process  are 
therefore  the  depth  of  cut,  feed  rate,  cutting  speed  and  the  width  of 
cut.   (There  are  also  other  parameters  such  as  the  rake  angle  of  the 
teeth  which  have  been  left  out  of  this  illustration  for  simplicity.) 
The  cutting  forces  generated  at  the  tooth  tip  may  be  resolved  into 
three  components,  viz.  the  radial  force  P  ,  axial  force  P  and  the 
tangential  force  P  . 
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OPTIMIZATION  METHODOLOGY 


The  mathematical  model  for  optimization  essentially  consists 
of  a  cost  objective:  function  subject  to  a  number  of  constraints  that 
arise  from  the  nature  of  the  physical  process  of  machining.  IjoJ 

The  total  unit  cost  of  the  milling  operation  consists  mainly 
of  three  parts:  machining  and  operating  costs,  tool  costs  and  non- 
productive costs.  The  machining  cost  per  unit  Is  the  cost  of  running 
the  machine,  cost  of  operator  and  overhead  computed  over  the  time  taken 
for  the  machining  of  a  single  piece.   The  tool  cost  is  made  up  of  tool 
setting  cost,  tool  regrlndlng  cost  and  tool  changing  cost  including 
cost  of  depreciation.   Finally,  the  non-productive  cost  includes  the 
cost  of  non-cutting  time  forced  by  the  nature  of  the  machining 
operation  and  is  equal  to  the  cost  of  operator  labor  and  overhead 
computed  over  the  non-cutting  time. 


17 
In  the  case  of  the  face  millinc  operation,  the  machining 

time  per  piece  t.  and  the  non-cuttinc  tice  t.  are  given  by 


dep      wid       T 

*-i        ^-     ''ij  "ij  "ij 


and 


fn,      n  ., 
dep     wid 


iJ  if  a  idle  return 

f-<      f-i         t; stroke  is  needed 

^-^     "J-^   .      ^rapid 


P  .  otherwise 


where  f  is  the  feed  per  tooth  (inches  per  tooth);  rapid  is  the  feed 

rate  of  the  work  table  for  rapid  traverse  in  the  return  stroke  (inches 

per  minute);  i,J  are  the  subscripts  used  to  denote  the  i'th  pass  in 

the  depth  direction  and  the  J'th  pass  in  the  width  direction;  L  is  the  length 

of  travel  of  the  cutter  relative  to  the  workpiece  in  each  longitudinal 

stroke  (inches);  N  is  the  number  of  Revolutions  per  minute  of  the 

cutter;  n  is  the  number  of  teeth  on  the  cutter  being  used;  n  .  ,  is  the 

number  of  passes  required  in  the  width  direction;  and  n^   is  the 

dep 

number  of  '^ccccs  in  the  dc^th  dirccticr.   Tt  tt""'   be  nctcd  th— t  L  ''s 

determined  by  the  longitudinal  dimension  of  the  surface  to  be  milled  and 

the  allowance  needed  for  the  cutter  approach;  n,   is  determined  by  the 

total  thickness  of  material  to  be  removed  and  the  depth  of  cut  given  in 

each  pass  in  the  deoth  direction;  and  n  . ,  is  determined  by  the  width 
*^  -  '      wid  '' 

dimension  of  the  milled  surface  and  the  width  of  cut  used  in  each  pass  in  the 
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width  direction.   The  total  cost  of  machining  can  be  determined  from 

a  knowledge  of  t-.t^  and  the  relevant  cost  rates  and  the  number  of 

tool  resettings  and  regrindings  required  for  the  machining  of  a 

single  piece.   It  should  be  mentioned  that  the  expected  number  of 

times  the  cutter  is  to  be  regrounded  and  reset  during  the  machining 

of  a  single  piece  is  given  by  E[t^/T]  where  T  is  a  stochastic  variable 

devoting  the  life  of  the  cutter  (i.e.,  the  time  it  can  be  used  for 

machining  before  failure  occurs),  and  E[    ]  stands  for  the  expectation 

operator.  Empirical  formulas  expressing  E(T)  as  a  function  of  the 

machining  parameters  are  available  in  literature  [10]  and  are 

generally  of  the  form  E(T)"  ^ —  where  C  and  Indices  x,y,z  are 

jX^yj^z        y 

constants  determined  by  the  tool  and  work  materials  being  used.  For 
the  present  purposes,  we  shall  approximate  the  expected  number  of  tool 
failures  (in  a  period  of  length  t^^)  by  t,/E(T). 

The  total  unit  cost  of  machining  can  be  written  as: 


+  (tj/I)   [(C3«  ^)t^^  +  (S^chg^chgl 


The  approximation  holds  very  well  for  the  case  where  tool  failures 
are  Poissori  distributed  and  t^  is  fairly  large  relative  to  E(T). 
Both  of  these  are  very  reasonable  assumptions  for  most  machining 
situations. 
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•where  C. .C-.C-  are  the  rates  for  direct  labor,  operating  overhead 

and  fixed  overhead  respectively,  t   ,t   ,t  .   are  respectively  the 

times  required  for  initial  tool  setting,  tool  regrinding  and  tool 

changing  (including  resetting)  and  C  ^,C  ,C  ,   are  the  wage  rates  for 

set  gr  cng 

tool  setting,  regrinding  and  changing. 

The  total  unit  cost  of  machining  should  be  optimized  with 
respect  to  the  machining  parameters  such  as  the  cutting  speeds,  feeds, 
depths  and  widths  of  cut.  However,  the  choice  of  these  machining 
parameters  is  restricted  by  a  number  of  constraints  that  result  from 
consideration  of  the  physical  process  of  milling  and  the  limitation 
of  the  available  facilities.  For  example,  the  spindle  speed  and  feed 
used  must  lie  within  the  range  of  maximum  speeds  and  feeds  available 
on  the  machine  tool;  the  cutting  forces  should  not  exceed  the  maximum 
forces  that  can  be  resisted  by  the  work  and  tool  holding  devices;  the 
cutting  power  required  should  be  less  than  the  net  power  available 
from  the  driving  motor  of  the  machine  and  so  on.  Certain  restrictions 
are  also  to  be  Imposed  to  Insure  chatterless  machining  and  to  limit 
the  probability  of  cutter  failure  during  the  machining  of  a  single 
surface  to  a  specified  maximum. 

The  major  constraints  for  the  case  of  face  milling  are  stated 

2 
below.   For  notatlonal  simplicity,  the  subscripts  l,j  (denoting  the 

sequential  number  of  the  pass  in  the  depth  and  vd.dth  directions)  are 


This  restriction  is  imposed  because  failure  of  the  cutter  during  the 
machining  of  a  single  workpiece  is  extremely  inconvenient  and  expensive 
in  terns  of  tine  loss  in  resetting  the  machine  and  cutter  in  the 
middle  of  a  pass. 

Once  again,  it  should  be  noted  that  the  constraints  may  look  somewhat 
different  for  other  machining  processes,  but  will  be  of  the  same 
general  forms. 
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omitted.  Except  when  explicitly  specified,  the  following  constraints 
apply  to  each  Individual  pass. 

Speed  Restrictions; 


where  N  .  and  N    are  the  minimum  and  maximum  available  spindle  speed: 
mln      max 

on  the  machine,  respectively.   There  are,  in  addition,  restrictions 
on  the  surface  speed  of  the  cutting  edge  relative  to  the  workplece. 
The  peripheral  speed  is  given  by  V»  (ttDN/12),  The  surface  speed 
constraints  can  then  be  written  as: 


12  V  .        12  V 
Si2.<  N  < ES2L 


i:D        ~  irD 


where  V    is  the  minimum  cutting  speed  to  avoid  the  formation  of 
built  up  edge  on  the  tool  tip.  V    is  the  maximum  cutting  speed  for 
the  prevention  of  excessive  heating  of  the  tool  tip,  and  D  is  the 
cutter  diameter  in  Inches. 
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Feed  Restrictions; 


F  .  <  fnN  <  F 

mln  -     ~  max 


where  F  .  and  F    are  the  minimum  and  maximum  feed  rates  available 
mln      max 

on  the  machine. 

Depth  of  Cut  Restrictions; 

d  <  H 

1   d  <  d 

-  max 

where  H  Is  the  total  depth  of  material  to  be  removed  after  correcting 

for  finishing  allowances  and  d    Is  the  length  of  the  cutting  edge 

max 

of  each  tooth  of  the  cutter. 

In  addition,  the  total  depth  of  material  removed  in  the  n. 
layers  should  equal  H.  i.e.. 


E   d.»  H 
i-1   ^ 
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Width  of   Cut  Restrictions; 


The  width  of  cut  w  Is  limited   to  a  icaxlmum  value  w . 

max 

For  most  face  milling  applications,  w  j^"  D  so  that 


w  <  D 


In  addition,  the  total  width  of  metal  removed  should  equal 
W,  the  width  of  the  workpiece.  Hence, 


r  „  -  W 
J-1  J 


Force  Restriction: 


The  maximum  resultant  cutting  force,  P   ,  at  any  tool  tip 

must  not  exceed  the  maximum  force  that  can  be  resisted  by  the  tool  tip, 

P   .  In  the  case  of  face  milling, 
max 


P   - 
res 


/P^T7< 


where  P  and  P  are  the  tangential  and  radial  cutting  forces  at  the 
tool  tip,  respectively.  The  radial  force  in  the  case  of  face  milling 
is  found  to  be  approximately  0.3  to  O.A  times  the  tangential  cutting 
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force  [4,12].   The  tangential  cutting  force  is  generally  computed  by 
.empirical  fomulas  expressing  it  as  a  function  of  each  of  the  machining 
parameters.  One  such  formula  knovm  as  the  extended  cutting  force  law 
la  given  by  Kronenberg  [10]  and  is  of  the  form 


P^-  C  Af  G^  N^ 
t   p  1 


%ifaere  k.    is  the  modified  area  of  the  chip  cross  section  and  is  given 
by  A."  df ,  G  is  the  slenderness  ratio  and  is  given  by  G=  d/f ,  and 
g,  h,  and  k  are  indices  which  depend  upon  the  particular  toolwork 
combination  being  used. 

If  P  -  C-P  ,  the  cutting  force  restriction  becomes. 


(^) 


x,  y,  z, 

C  d  -^f  -^N  -^  <  P 

p  I 


where  x-,  y.  and  z.  are  indices  determined  by  the  values  of  g,  h,  and  k. 


Horsepower  Restriction;  -,        -        _ 

The  total  power  available  at  the  cutter  must  be  enough  to  resist 
the  cutting  forces  at  the  speeds  being  used.   This  requires  that 


P  V 

t-tot  ^   .     , 

33000        -  ^^^   Vach  Vjt 
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where  HP  is  the  horsepower  of  the  driving  motor,  P-^_^  ^  is  the  total 

of  the  tangential  force  components  with  all  the  teeth  in  contact  with 

the  workpiece,  n   >,  is  the  efficiency  of  power  transmission  in  the 

machine,  n  ^  is  the  efficiency  of  the  driving  motor,  and  V  is  the 
mot 

peripheral  velocity  of  cutting.  P^_,^  f.   assumes  its  maximum  value 
when  all  the  teeth  of  the  cutter  are  in  contact  with  the  work  surface, 
and,  in  this  case,  it  is  given  by 


t-tot  2  ""t 

Substituting  this  and  the  expression  for  peripheral  speed  V,  the 
power  constraint  becomes- 


C     _„   X-  y-  z.+l 
-E^C^nd^fHl^  )<HPn. 


33,000  ^24  ""         /  ^  •-  -jjj^^^  ..jjj^^ 

Restrictions  to  Avoid  Excessive  Vibrations  and  Chatter; 

Chatter  is  a  function  of  both  the  cutting  conditions  and  the 
machine  structure;  thus  it  is  not  easy  to  form  constraints  that  will 
positively  insure  chatterless  operation.   Based  upon  experimental  and 
theoretical  research  of  Tobias  [22],  Tobias  and  Sadek  [23],  Lemon  and 
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Long  [11],  and  others,  the  following  constraint" is  applied 


=^-v 


where  n  is  the  number  of  edges  cutting  simultaneously,  X  is  the 
c 

directional  static  stiffness  of  the  machine  frame  (Ibf/in.),  k^.  is 
the  slope  of  cutting  speed-cutting  force  curve,  and  k^  is  a  constant 
parameter  indicating  the  safe  region  with  respect  to  vibrations. 
(A  typical  value  of  k  Is  0.2).  The  maximum  value  of  n^  Is  n/2  and 

^t        ^1  ^1  ^1'^ 
k  -  -f--  C  x-d  ^f  -^N  ^ 
c    d     p  1 


so  that  the  chatter  constraint  becomes 


/i/i/r^    2k^x 
p 


Other  researchers  like  Barash  [2]  suggest  consideration  of 
the  chatter  problem  by  a  constraint  on  the  chip  slenderness  ratio.  In 
this  work,  the  constraint  used  is 

d/f  <  50 

Hopefully,  these  two  constraints  used  together  will  keep  the  machining 
parameters  in  the  region  permissible  with  respect  to  chatter. 
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Restriction  of  Tool  Life; 

It  is  highly  desirable  that  a  freshly  reground  cutter 
last  for  at  least  the  time  required  to  machine  a  single  workpiece. 
Failure  of  the  cutting  tool  in  the  middle  of  machining  is  extremely 
Inconvenient  and  expensive  in  terms  of  time  lost  in  resetting  the 
machine  and  cutter  in  the  middle  of  a  pass.  However,  since  tool 
failure  is  inherently  probabilistic  in  nature,  this  constraint  can 
never  be  satisfied  with  absolute  certainty.  Hence,  it  is  converted  to 
a  * chance  constraint'  of  the  form 

PrfT  >  tj^}  2  a. 

where  a  is  chosen  depending  on  the  degree  of  assurance  with  which  the 
tool  life  constraint  is  desired  to  be  satisfied.   It  is  obviously  in- 
convenient to  handle  the  tool  life  constraint  in  its  probabilistic  form 
directly.   If  the  probability  distribution  of  the  tool  life  is  known, 
the  constraint  can  be  transformed  into  a  more  convenient  form  by 
using  the  properties  of  the  distribution.  The  foregoing  inequality  can 
be  rewritten  as 


Pr  {- ^^^^  >  — }  >  a 
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vhere  E(T)  and  a_  are  the  mean  and  standard  deviation  of  the  tool  life 
respectively.   Let  K  be  < 
which  T  belongs  such  that 


respectively.   Let  K  be  a  constant  defined  for  the  distribution  to 


Pr{e  >  K  }-  a 
a 


where  e  is  the  standardized  variate  of  the  distribution.   Using  this 
notation, 


Pr{i^-Effl->K 


.>■» 


life  constraint  can  thus  be  expressed  as 


t,  -  E(T) 

o^ -^« 

T 


Substituting  for  t-  and  transforming, 


^dep  '^wid      L  C' 

TV  ^j  y 


<  K 


i-1  J-1    f^jN^j  n,j   -  ,x,y^z  i  \'T 

Very  little  experimental  work  has  been  done  to  establish  the  nature  of 
the  probability  distribution  of  tool  life.  It  is  known,  however  [21], 
that  two  different  t>T3es  of  tool  failure  are  prominent  in  the  case  of  . 
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carbide  tools,  failure  due  to  excess  flank  wear,  and  failure  due  to 
chipping.   It  can  be  expected  that  the  probability  distributions  of 
life  of  the  tool  would  be  different  for  the  two  different  modes  of 
failure.  The  Poisson  distribution  has  been  suggested  for  tool 
failures  by  chipping  [1]  and  some  form  of  modified  normal  distribution 
for  the  tool  failure  by  flank  wear.  The  actual  distribution  of  tool 
life  would,  therefore,  be  the  composite  distribution  obtained  by 
taking  the  weighted  mean  of  the  distributions  in  the  case  of  the  two 
different  modes  of  failure,  welghtages  being  assigned  on  the  basis  of 
relative  dominance  of  each  mode  of  failure.  Moreover,  the  failure  of 
a  milling  cutter  is  governed  by  the  failure  of  any  one  of  its  teeth 
and  so  the  appropriate  distribution  is  the  distribution  of  the  minimum 
of  the  lives  of  the  individual  teeth  of  the  cutter.  Therefore,  the 
probability  distribution  of  the  life  of  a  cutter  has  been  derived  from 
that  of  an  Individual  tooth  which  in  turn  is  obtained  from  the 
distributions  of  failures  by  chipping  and  failures  by  flank  wear.  With 
the  distribution  of  T  known,  one  can  easily  read  off  the  value  of  K 


for  any  chosen  value  of  a. 
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FINAL  FORM  OF  THE  OPTIMIZATION  MODEL 


The  optimization  problem  as  stated  above  is  quite  a  formidable 
one,  one  with  non-linear  objective  function  and  many  non-linear  constraints. 
By  a  detailed  consideration  of  the  nature  of  the  objective  function,  nature 
of  the  constraints  and  the  physical  process  of  milling  the  reasonable 
modifications  and  simplifications  of  the  above  mathematical  mod^l  have 
been  made.    It  is  also  found  convenient  to  define  the  transformed  variables 

V  1°8  (^/^mln) 
and  X3-  log  (N/N^^^). 

Expressed  in  terms  of  those  transformed  variables,  the  above 
mathematical  model  can  be  shown  to  reduce  to  a  convex  objective  function 
with  a  set  of  linear  constraints.   It  reduces  to  the  form 


„,  ,  ,   ,   *=10^1lV^12'^2-'^13^3  ^  «=20-^2lV*=22'^2+'^23'^3 
nlnlnlze  Z"  e  +  e 


subject  to  a  number  of  constraints 


ax  +  bj^  <  0    1=  1 1  and  j=  1,2,3. 

X.  >  0 


where  the  b.  s  and  c..  s  are  functions  of  the  constants  and  parameters  of 
the  model.  . 
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The  technique  of  optinlzatlon  employed  is  a  version  of  the 
gradient  projection  method  very  similar  to  that  suggested  by  Rosen  [17], 
This  method  belongs  to  the  general  class  of  'hill-climbing'  methods 
for  optimization  of  non-linear  functions.   It  always  attempts  to  travel 
In  the  locally  steepest  direction  of  ascent  (or  descent)  while  not 
allowing  the  point  to  leave  the  constrained  set.   The  optimum  size  of 
the  step  is  then  determined.   Once  the  optimum  Is  found  on  the  line 
of  steepest  ascent,  a  new  gradient  evaluation  is  carried  out  and  the 
procedure  repeats  itself  until  a  point  sufficiently  close  to  the  global 
optimum  is  reached. 

Notice  that  the  objective  function  we  are  encountered  with  is 
very  smooth,  convex  and  'well-behaved'  and  satisfies  all. the  regularity 
conditions  needed  for  convergence.   It  is-  therefore  reasonable  to 
expect  the  procedure  to  converge  in  a  relatively  small  number  of 
iterations.  Moreover,  the  objective  function  is  composed  of  the  sum 
of  two  exponential  terms  which  makes  it  extremely  easy  to  analytically 
compute  the  gradient  vector  at  any  stage  of  the  iterative  procedure. 
It  is  precisely  these  features  that  make  the  gradient  methods  so 
convenient  for  this  class  of  problems. 

The  optimization  procedure  is  carried  out  by  means  of  the 
computer  program  OPTIMIZE  which  is  part  of  the  software  that  comprises 
the  planning  and  optimization  system. 


This  is  by  no  means  an  accident;  nor  is  it  peculiar  to  the  face  milling 
process.   It  can  be  shown  that  the  model  will  reduce  to  a  similar  form 
for  most  machining  processes. 


AN  ILLlI3TRATT0:i  OF 


USE  OF  THE  SYSTE:-! 
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The  use  of  the  suggested  system  is  illustrated       !. 
b«low  for  the  case  of  a  very  simple  milled  part  shown  in 
Figure  3.  It  is  asstimed  that  the  raw  material  selection  phase  of 
the  planning  has  been  conpleted  ajid  the  system  has  chosen  the  rav 
material  shovn  in  Figure  h.     The  dimensions,  tolerances  and  surface 
finish  specifications  of  the  finished  component'  are  given  by  the 
Component  Data  iMatrix  (CDM)  which  is  reproduced  in  Table  1.   The 
machining  of  the  part  has  been  divided  into  five  stages  of  milling, 
each  stage  corresponding  to  milling  of  a  plain  surface.  The  five 
rows  of  the  CDM  thus  represent  the  dimensions  and  specifications 
of  these  five  plain  surfaces  as  they  should  appear  after  each 
stage  of  milling. 

■  TABLE  I 

•      '  •'.    COMPONENT  DATA  I-L\TRIX 


LEIs'GTHS  (INCHES) 

TOLERANCES  (THOUa) 

SURFACE  ■      1 
FINISH  GRADES 

NO 

LX 

LY 

LZ 

TOLX 

TOLY 

TOLZ 

SFX 

SPY 

SFZ 

1  • 

lU.OO 

5.00 

8.75 

999.00 

999.00 

50.00 

0 

0 

h 

2 

8.75 

5.00 

12.50 

999.00 

999.00 

30.00 

0 

0 

5 

3 

12.50 

8.75 

I4.OO 

999.00 

999.00 

10.00 

0 

0 

8- 

k 

12.50 

I4.00 

2.50 

999.00 

999.00 

10.00 

0 

0 

5 

5 

I.. 57 

6.00 

3.00 

999.00 

999.00 

5.00 

0 

0 

•  10 

Plgurc  3.  .  Finished  Workpiece 


Figure  A.  Rav  Workpiece 
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The.  dincnclonc  of  the  raw  natarial  for  each  stage  of 
milling  are  given  by  the  Raw  Material  CharacterisCics  Matri::  (KMCM) 
shown  In  Table  II. 

•TABLE  II  .     . 

RAW  MATERIAL  CHARACTERISTICS  MTRIX        ''        _.  '. 


WORK  SURFACE  NO. 

LX 

LY 

LZ 

1 

lA 

5 

10 

2 

8.75 

5 

lA  ■ 

3 

12.5 

8.75 

5 

4 

12.5 

A 

A 

5 

A. 75 

6 

A 

It  is  assumed  In  this  example  that  only  two  machines  and  two  cutters 
are  available  and  their  specifications  are  glyen  by  the  Machine 
Characteristics  Matrix  (MC7I)  and  the  Tool  Characteristics  Matrix  (TCM) 
sho^m  in  Table  III.   In  the  Machine  Characteristics  Matrix,  column 
headings  ISIIi;  and  'HL'iX  stnr.d  for  the  r.ininum  ar.d  Tnnxinun  r-pindlp.  Knec>ds 
available  on  the  machines,  FMIN  and  FMAX  are  the  ninimua  and  naxirsuE 
feed  rates  available,  HP  is  the  horsepo;7er  of  the  driving  motor,  EFFf  1 
and  EFPy  2  are  the  efficiencies  of  the  driving  motor  and  the  trans- 
alsslon  of  the  machine  respectively,   LAMBDA  is  the  stiffness 
of  the  machine  frame,  CI,  C2,  C3  are  the  different  wage  and  overhead 
rates  and  TSET  is  the  time  required  for  setting  the  machines. 


In  the  case  of  the  Tool  Characteristics  Matrix,  the  coliu?.ns  read 
respectively,  cutter  number,  nur.ber  of  teeth,  diameter  of  the  cutter,. 
length  of  the  cutting  edge,  naxir.iun  force  that  can  be  resisted  by 
the  tooth  tip,  tine  for  regrinding  and  tin-.e  for  changing  the  cutter. 
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TABLE  III 

f4AC}{INE  CHAR.'VCTKRISTICS  MATRIX 


NO. 


NMIN 


NMAX 


FMIN 


FMAX 


HP 


EFFYl 


JJ'FY2 


LA>II5DA 


CI        C2        C3      TSI 


30.0 
AO.O 


800.0 
500.0 


0.2 
0.3 


5.0 
3.0 


20.0 
15.0 


0.7 
0.7 


1.0 
0.9 


210000.0 
1]0000.0 


0'.2 
0.3 


0.2 
0.2 


1.1  20. 


0.1 


AO. 


TOOL  CHAR/\CTERISTICS  Il.^TRIX 


NO. 

TEETII 

DIA 

LEDGE 

PMAX 

TGR 

TCHG 

1 
2 

■8 
12 

4.0 
6.0 

0.4 
0.2 

LOOO.O 
2000.0 

190.0 
135.0 

75.0 
60.0 

The  program  MODIFY^  accepts  the  CDM,  RMCM,  MCM  and  TCM  as  data, 
dctcnnir.cs  the  v:iric'jc  Tnachining  requirements  on  each   vovk  surface 
and  generates  the  Machine  Requirements  Matrix  (MRM)  described  earlier 
In  the  article  as  well  as  the  matrices  CDMRG,  CDMFC,  CDMRC  and  RCSIZE. 
These  matrices  are  reproduced  in  Table  IV. 


All  the  progra.Tis  named  here  are  parts  of  the  software  that 


make  up  the  system. 


TABLE   IV 

machii:e'  reouipz:z:jt 


^ix 


SURFACE  i;a 

RC 

FC 

RG 

EG 

0 

0 

0 

1 

0 

0 

1 

1 

0 

1 

1  • 

0 

1 

1 

1 

RDf^RC 


14. 

5 

8.8 

8.75 

5 

12.57 

12.5 

8.75 

4.08 

12.5 

A 

2.58 

A. 75 

6 

3.085 

CDIvIFC 


1 

14 

5 

8.75 

2 

8.75 

5 

12.52 

3 

12.5 

8.75 

4.03 

4- 

12.5 

4 

2.53 

5 

4.75 

6 

3.035 

CDI-IRG 


1 

14 

5 

8.75 

2 

8.75 

5 

12.5 

3 

12.5 

8.75 

4.01 

4 

12.5* 

4 

2.51 

5 

4.75 

6 

'.  3.015 

CDMFG 


1 

14 

5 

8.75 

2 

8.75 

5 

12.5 

3 

12.5 

8.75 

•  4. 

4 

12.5 

A 

2.5 

5 

4.75  . 

6 

3.005   1 

ROUGH  CUT  SIZE  MS.TRIX 


1 

14 

5 

1.2 

2 

8.75 

5 

.1.43 

3 

12.5 

8.75 

0.92 

4 

12.5 

4 

1.42 

5 

4.75 

6 

0.915 
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STAGE  1 
SCALE  10:1 


8,75 


STAGE  2 
SCALE  10:1 


cp 


Ui 


■12.5- 


8.75 


STAGE  3 
SCALE  10:1 


STALE  i+ 
SCALE  10:1 


^xO^^^.-^   7^   J   STAGE  5 
f^/r^;,,,,,.-  ^"^     <  J(U    SCALE  10 

^.75 


Figure  5.  Optimum  Sequence  of  metal  removal  for  the  sample  part. 
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The  progran  CALCULATE  then  uses  the "data  generated  by  MODIFY 
to  calculate  the  values  of  the  various  coefficients  and  constants 
In  the  constraint  matrix  and  the  objective  function.   The 
program  CHOOSE  considers  all  the  available  machine-cutter 
combinations  (4  in  this  case)  one  after  the  other  and  uses  the 
programs  MODIFY,  CALCULATE  and  OPTIMIZE  to  determine  tlie  optimum 
conditions  for  each  combination.   It  then  chooses  the  best 
combination  of  machine  and  cutter  and  the  corresponding  optimized 
machining  parameters.   These  results^,  along  with  the  nuriber 
of  passes  required  in  the  depth  and  width  directions  and  estimates 
of  total  machining  time  and  machining  cost  per  piece  are  printed 
out  by  the  system.   For  purpose  of  making  comparisons  and  any 
necessary  compromises,  the  optimized  costs  of  machining  with  all 
other  possible  machine-cutter  combinations  also  are  printed  out 
*  as  supplementary  data.   The  reconunendations  of  the  machines, 
cutters  and  machining  parameters  for  each  of  the  five  work  surfaces 
of  the  present  example  are  presented  in  Table  V. 


TABLE  V 

RECo:rinr:D  MAC!iT:;ns,  cutters  ,v:n  yu.\cHTNiNr.  parameters 
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WORK 

M\CHINE 

CUTTER 

CUTTER 

NUMBER 

RPM  OF 

FEED 

DEPTH 

SURFACE 

I.'UMEER 

lOJMUER 

DI.\:iETER 

OF  TEETH 

CUTTER 

ao-3 

OF  CUT 

WUllBCR 

(INCHES) 

IN   CUTTER 

INCHES 
PER  TOOTlv 

(INCHES) 

1 

1 

1» 

8 

71 

5.7 

O.2I4O 

2 

2 

6 

12 

80 

U.O 

0.176 

3 

2 

6 

12 

80 

U.O 

0.18U 

4 

1 

k 

8 

71 

5.7 

o.2eu 

5 

2 

6 

12 

80 

1*.0 

0.183 

WIDTH 
OF  CUT 
(INCHES) 

"dep 

"wid 

ESTIMATED 
]4ACHIIIi:iG 
TIME  PER 
PIECE    (MIN.) 

ESTIMATED 
MACHINING   COST 
PER  PIECE 

2,500 
5.000 
I..  375 

U.ooo 

6.000 

2 

1 
2 

1 
1 

5 
8 

5 
5 
5 

kh,9U 

22.02 

36.8 

23.8 

16.7 

36. 6U. 

22. 9^ 
31.36 
22. 5^* 
11.01 

CONCLUDING  REMARKS 


The  system  concept  presented  herein  is  general  enough  to  be 
applicable  to  a  variety  of  machining  processes  and  diverse 
machining  environments.      Of  course,  the  detailed  formulation  of 
the  planning  and  optimization  models  must  be  developed  keeping 
in  mind  the  features  of  the  specific  machining  situation.  For  e:^anple, 
the  minimum  cost  criterion  has  been  used  in  the  illustration  given 
here,  but  this  may  or  may  not  be  the  most  appropriate  criterion  in 
a  specific  machining  problem.   Other  authors  have  suggested  various 
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other  criteria  such  as  maximum  production  rate  [1]  maximum  profit  [13], 
minimum  time  [3],  etc.   Similarly,  while  working  on  automatics 
employing  several  tools  simultaneously,  fixed  tool  servicing  schedules 
might  be  preferable  to  changing  the  individual  tools  when  failure  is 
expected. 

The  authors  believe  that  systems  for  automated  manufacturing 
planning  such  as  the  one  suggested  in  this  article  are  likely  to  find 
wide  practical  utility  in  the  near  future.  Manufacturing  systems  are 
typically  complex  and  extremely  djmamic.  Development  of  such  ' self- 
optlmizing '  units  should  prove  to  be  a  major  technological  advance  for 
the  manufacturing  industry.  , 
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